Two-dimensional Ising model with self-dual biaxially correlated disorder 
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We consider the Ising model on the square lattice with biaxially correlated random ferromagnetic 
couplings, the critical point of which is fixed by self-duality. The disorder, which has a correlator 
G(r) ~ r~ , represents a relevant perturbation according to the extended Harris criterion. Critical 
properties of the system are studied by large scale Monte Carlo simulations. The correlation length 
critical exponent v = 2.005(5) corresponds to that expected in a system with isotropic correlated 
long-range disorder, whereas the scaling dimension of the magnetization density x m — f3/v = 
0.1294(7) is somewhat larger than in the pure system. Conformal properties of the magnetization 
and energy density profiles are also examined numerically. 



I. INTRODUCTION 

The presence of quenched disorder of different types is 
an inevitable feature of real materials. The effect of ran- 
domness on the physical properties of systems is often 
dramatic, in particular in the vicinity of singular points, 
such as for classical phase transitions at finite temper- 
ature or for quantum phase transitions at T = 0. In 
a random system the set of critical exponents associ- 
ated with a phase transition can be completely differ- 
ent from that in the nonrandom system, which means 
that disorder can be a relevant perturbation in the sense 
of the Harris criterion^ In most of the studies of dis- 
ordered systems the parameters (couplings, longitudinal 
or transverse fields, etc.) are expected to be indepen- 
dent and identically distributed random variables among 
which there is no correlation. 

There are, however, systems in which the random vari- 
ables exhibit some kind of correlations, the effect of which 
can even change the random universality class of the 
phase transition. Examples of systems with strongly cor- 
related disorder are random quantum magnets at T = 0, 
in which the disorder is strictly correlated in the tem- 
poral direction. Indeed the critical properties of a such 
type of random systems are quite unusual, for a review 
see Ref. 0. In the classical counterpart of these quan- 
tum systems, the disorder is strictly correlated in one di- 
rection, which leads to models introduced and partially 
exactly solved by McCoy and Wui^ In a further general- 
ization, the disorder can be strictly correlated in a d' < d 
dimensional subspace, where d is the dimension of the 
system. 4 i 5 i 6 i 7 i 8 

Correlations between random variables can also be 
isotropic. One example is presented in Ref. |9| in which 
the experimental observations are explained^ with the 
assumption of the presence of randomly oriented dis- 
location lines in the sample. Also in some random 
quantum magnets** and in the non-Fermi liquid be- 
havior of /-electron compounds 1 — the disorder is ex- 
pected to be isotropically correlated due to long-ranged 



Ruderman-Kittel-Kasuya-Yoshida (RKKY) interactions. 
In a coarse-grained picture, the disorder in some coupling 
strength AK(r) = K(r) — K(r) is characterized by its 
correlator, G(r — r') = AK(r)AK(r'), which is expected 
to be isotropic, G(r) = G(|r|), and to decay algebraically 
for large arguments: G(r) ~ r~ a . For example, randomly 
oriented straight dislocation lines are associated with a 
correlator exponent a = d — 1. 

The influence of Gaussian disorder with power-law cor- 
relations on the critical behavior of the m-component 
vector spin model has been studied in Ref. by field- 
theoretical methods using a double expansion in e = 4 — d 
and 8 = 4 — a. Extending the Harris criterion, correlated 
disorder with a > d is found to be effectively short range 
and thus relevant only when the pure system has a pos- 
itive specific heat exponent i* When a < d and m = 1 
(Ising model) the correlated disorder is shown to be a 
relevant perturbation at the short-range-disorder fixed 
point, provided the correlator decays sufficiently slowly, 
i.e., when 



a < 



fshort 



(1) 



where z4hort is the correlation length exponent of the sys- 
tem with short-range disorder. At the stable long-range- 
disorder fixed point the correlation length critical expo- 
nent is found to be given by 



= 2/a, 2/a > u, 



short - 



(2) 



This result is argued to be exact and thus the extended 
Harris criterion in Eq. is marginal when the long- 
range nature of the correlations is relevant. On the other 
hand the magnetic exponents, such = (3/v, depend 

on both a and d. 

Recently^ the validity of the scaling relation J2J has 
been questioned on the basis of a direct renormalization 
analysis of the scaling functions for the three-dimensional 
(3D) system in the two-loop approximation, using a 
Pade-Borel summation technique to evaluate the critical 
exponents for different values of a. 
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Numerical and field-theoretical calculations on dif- 
ferent types of systems with correlated disorder 
(percolation^ 3D Ising model^ polymers^ quantum 
Ising model^i etc.) seem to be in overall agreement with 
the prediction of Eq. J5J. 

Here we should make, however, two remarks. First, in 
the case of the Ising model with randomly distributed dis- 
location lines, oriented along the lattice axesji^ for which 
rotational symmetry is broken, the critical behavior was 
found to be the same as for the model with Gaussian cor- 
related disorder when a = 2. This shows that the details 
of the model of correlated disorder are seemingly irrele- 
vant, the critical behavior being essentially governed by 
the decay of the disorder correlations. Our second re- 
mark concerns the fact that, in systems with correlated 
disorder, the location of the critical point is generally 
not known exactly, which introduces a limitation on the 
accuracy of the numerical calculations. 

In the present paper, our aim is to increase the nu- 
merical accuracy of the calculation in a statistical me- 
chanical system with correlated disorder and address also 
such questions (c.f. critical point density profiles, valid- 
ity of conformal invariance, etc.) which have not been 
considered previously. For this purpose, we consider 
the two-dimensional (2D) square lattice Ising model and 
introduce a type of biaxially correlated disorder which 
preserves self-duality, thus the critical point is exactly 
known. Our model can be considered as a 2D self-dual 
version of the random dislocation model of Ballesteros 
and Parisi^ In the vicinity of the critical point the 
properties of the system are explored at the scale of the 
(diverging) correlation length. A correlator exponent, 
o = l, can be deduced from the integral of the disorder 
correlator inside a large square of linear size, L. Thus 
at a coarse-grained level our model can be compared to 
a model with Gaussian correlated disorder for the same 
value of a. 

For the 2D Ising model uncorrelated disorder coupled 
to the local energy density (such as dilution or ran- 
dom ferromagnetic couplings) is a marginally irrelevant 
perturbation^*^ hence critical singularities of the pure 
system are supplemented by logarithmic corrections. In 
particular the correlation length exponent J4hort keeps 
its unperturbed value vq = 1 . As a consequence from 
Eq. (JJJ the biaxially correlated disorder of our model is 
expected to be relevant. In order to explore the criti- 
cal properties of the system and to check the validity of 
Eq. J5J we perform large scale Monte Carlo simulations. 
In particular the critical exponents are deduced from the 
finite-size scaling behavior of various magnetic and ther- 
mal quantities at the critical temperature, which is ex- 
actly known. We also determine the critical profiles of 
the magnetization and energy densities in order to check 
whether the forms obtained for the profiles through con- 
formal methods remain valid in the random system after 
averaging over the disorder. 

The structure of the paper is the following. The model, 
its self-duality, and the relevance- irrelevance criterion are 



described in Sec.[n] Finite-size scaling calculations of the 
critical exponents are presented in Sec. lIIII and the critical 
profiles are studied by conformal methods in Sec. II VI Our 
results are discussed in Sec. 

II. MODEL OF SELF-DUAL BIAXIALLY 
CORRELATED DISORDER 

A. Random-bond Ising models 

We study the spin 1/2 2D Ising model on a square 
lattice with Hamiltonian 

L L 

-PH = E E ( K fj^i+hj + Kfj^id+i) , (3) 
i=i j=i 

where (3 = l/fc^T. The Ising spin Uij = ±1 is associated 
with the site located at the intersection between 

column i and line j. Kfj and Kfj denote the random 
couplings between <7jj and its first-neighbors Cj+ij in the 
horizontal direction and <Ji,j+i in the vertical direction, 
respectively. In the following we use the parametrization 

e 2K h-l = 2 1 / 2+u h, (4) 

with t — x,y. Under duality transformation^ the hori- 
zontal and vertical directions are exchanged and we ob- 
tain a simple relation between the ufj and uf^ variables 



where the superscript, *, is used to denote dual variables. 

With different distributions of ujj and thus K[j differ- 
ent models are defined, which we list below. 

1. Nonrandom model 

The nonrandom model with uf, = u x and uf„- = u y has 
its critical point at u x = —u y , which follows from self 
duality in Eq. (JSJ). For the isotropic model the critical 
point is at u x = u v = 0. The critical exponents are 20 
v = 1 and x m — f3/i> = 1/8. 

2. Random model with uncorrelated disorder 

In this case the ujj parameters are independent and 
identically distributed (iid) random variables, which are 
taken from a distribution P(u). At the critical point 
of the model the distribution is symmetric, 21 P(u) = 
P(—u), which is also a consequence of self-duality, see 
Eq. The strength of disorder A is measured by 

A 2 = u 2 . In a coarse-grained picture cells, of linear size 
L, are defined with cell variables, v! = L~ d ^ Uij, where 



the sum runs over L sites, and d = 2 is the dimension 
of the system. In the coarse-grained description we have 
A' = AL~ d / 2 . To decide about relevance-irrelevance of 
a perturbation according to Harrisi one should consider 
the ratio A'/t at a length-scale, L = £ ~ t^", where t 
measures the distance from the critical point. For the 2D 
Ising model this analysis predicts a marginal effect of the 
disorder, actually a marginally irrelevant one. As a con- 
sequence critical singularities of the nonrandom model 
are supplemented by logarithmic corrections^*^ 

3. Random model with isotropically correlated disorder 

In this case the ujj parameters are random numbers 
which are correlated at different sites and we have for 
large spatial separation 

mjWJTr. [(i_i') 2 + (j-j') 2 p /2 . (6) 

The strength of disorder in the coarse-grained picture 
reads asi& A' c — AcL - ", with u = min(d, a)/2, thus 
long-range correlation in the disorder can be relevant for 
a < d. Indeed for the 2D Ising model with a < 2 this type 
of perturbation is relevant, see Eq. Q and the correla- 
tion length critical exponent is conjectured to be exactly 
given by Eq. (J2J. For isotropic correlated disorder the 
location of the critical point is not known exactly which 
makes it difficult to analyze numerical results in the crit- 
ical regime. 

4- McCoy-Wu model 

In the McCoy- Wu model^ the disorder is strictly cor- 
related in columns (or lines) 

u& = 6?(=aj), 4=af(=&p, (7) 

where the parameters 6f , a\ (aj, 6p are iid random vari- 
ables. Now in the coarse-grained description with cells 
of size L, the disorder strength behaves as A' MW = 
AmwL^ 1 / 2 . This type of disorder represents a relevant 
perturbation. According to exact and conjecturedly ex- 
act results, this model displays a strongly anisotropic 
scaling behavior, 22 so that ln£j_ ~ y^, where £ 
is the correlation length in the nontranslationally in- 
variant (translationally invariant) direction. The critical 
exponents^ are given by: v = 2 and x m = (3 — v5)/4. 

B. Biaxially correlated disorder 

Here we introduce a type of correlated disorder which 

is 

• Symmetric, with respect to the interchange of the 
x and y axis, 
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FIG. 1: Parameters entering in the definition of the correlated 
disorder. 

• Has a duality symmetry and 

• Represents a 2D version of the random dislocation 
models 

To be concrete we use the parametrization 

>>: ■ <<; • <- = <+&j, (8) 

(see Fig.^), which is the local average (for each bond) of 
the couplings in two McCoy- Wu models where the disor- 
der is along the x and the y axis, respectively, see Eq. (J7J . 
Consequently the strength of disorder is rescaled as in the 
McCoy- Wu model, A' b = A 6 L _1 / 2 , but the system has 
isotropic scaling, as far as the asymptotic behavior of the 
average correlations is considered. This type of pertur- 
bation is relevant (see Sec lll Cf) and has the same scaling 
exponent as the isotropically correlated random pertur- 
bation with the correlation parameter a = 1 in Eq. 10. 
We expect that the critical behavior is dictated by the 
averaged decay of the disorder correlations. As for the 3D 
random dislocation modeli& where the dislocation lines 
are located along the lattice axes, the breaking of ro- 
tational symmetry (before averaging over the quenched 
disorder) should not change the universality class. The 
present model, which has the same type of disorder corre- 
lations as the Ballesteros-Parisi dislocation model, should 
also display the same critical behavior as the model with 
isotropically correlated disorder discussed in Ref. IT3I 

An important advantage of the parametrization used 
in Eq. © is that the couplings obey the duality relation 
in Eq. J5J and thus the system is self-dual for an appro- 
priate choice of the coupling parameters a\ , , bf , and 
tij. Self-duality can be realized in different ways, leading 
to a critical system in which both lattice directions are 
statistically equivalent or not. Here we study the geome- 
try of disorder with the highest symmetry, i.e., with the 
same symmetric and uniform probability density for all 
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the random variables 



havior of the ratio 



P(s) = P(- fl ) 



1 

2A 



-A < s < A 
Isl > A 



a v a x b x 



I?. 



(9) 



The evolution of the average value K av and the standard 
deviation AK of the coupling as a function of the half- 
width A of the probability density is shown in Fig. [21 
The standard deviation, with the expansion 

AK = .165764A + .000085A 3 + .000010A 5 + 0(A 7 ) , 

(10) 

is almost linear in A. Evidently the relative strength of 
the disorder AK / K av is monotonously increasing with 
A. 



C. Relevance-irrelevance criterion 

The stability of the fixed point governing the critical 
behavior of the pure 2D system in the presence of corre- 
lated disorder can be analyzed as by Harrisii 3 - At the 
length scale L, the sum of the deviations from the average 
coupling grows typically like AL 3 / 2 since it contains 2L 
independent random variables with a vanishing average, 
each variable being of order AL due to the correlations 
along the L lines and the L columns. Dividing by the 
number of couplings one defines a typical deviation from 
the average temperature as 



r3/2 

5t(L) ~ A^— - ~ AL" 1 / 2 . 
L^ 



ill) 



The relevance of the perturbation depends on the be- 




10 12 14 



FIG. 2: Variation with A of the average value K av and the 
standard deviation AK of the couplings. The standard devi- 
ation is almost linear in A. 



t 



Co oc t 



(12) 



when one approaches the critical point (t — » 0). Here 
t = \T — T c \/T c is the reduced temperature, £o the cor- 
relation length, and vq the correlation length exponent 
of the unperturbed system. The perturbation is relevant 
when this ratio diverges at the critical point, i.e., when 
vq < 2. Thus the biaxial correlated disorder is a relevant 
perturbation in the case of the 2D Ising model for which 
v a = 1. 

Alternatively the same conclusion can be drawn from 
the value of the scaling dimension of A following from 
Eq. I|ll|) . Under a change of the length scale by a fac- 
tor b, St, which is a thermal perturbation, has a scaling 
dimension yto = l/v<j and thus transforms as 



St' = A'L'- 1 ' 2 = b^St = b^^AL- 1 ' 2 . 
With L' = L/b one obtains 

A' = b y& A, y A = 



2v Q 



(13) 



(14) 



so that, in agreement with the previous result, one finds 
that the strength of the disorder increases under rescaling 
when vq < 2. 

Here we remind that the biaxial correlated disorder, 
which is considered in this paper, can be represented 
by an effective correlator exponent a = 1 Indeed, the 
stability limit of the fixed point for Gaussian correlated 
disorder in Eq. Q with a = 1 is the same as the one ob- 
tained here for biaxial correlations. As mentioned in the 
Introduction, this result follows from the comparison of 
the integral of the disorder correlator at a length scale L 
for both models. The possible universality of the critical 
properties of the system with respect to the form of cor- 
related disorder, in particular the validity of the result 
predicted for the correlation length critical exponent in 
Eq. J2| will be studied in the following sections. 



III. FINITE-SIZE SCALING AT THE CRITICAL 
POINT 

A. Monte Carlo technique 

We study the finite-size scaling behavior of different 
magnetic and thermal quantities at the critical point. We 
work on a square-shaped system with size L and we use 
periodic boundary conditions in both directions. 

In order to limit the effects of the critical slowing down, 
the Monte Carlo simulations are performed using the 
Wolff cluster algorithm. 23 . In order to estimate the equi- 
libration time r e measured in cluster flips units, we com- 
pared the evolution of the energy and the magnetization 
towards equilibrium starting either from a random ini- 
tial state or from an ordered initial state. We have also 



5 



studied the autocorrelation time for the energy and the 
magnetization at equilibrium, in order to evaluate the 
statistical errors on the time averages. 

The time averages are calculated during a time r = 10 5 
to 2.5 x 10 5 Monte Carlo steps (MCS) after a waiting time 
of the order of 10r o (10 3 to 5 x 10 4 MCS). The disorder 
averages are taken over n s « 10 4 samples. These values 
were chosen in order to obtain a precision close to 1% for 
the different quantities studied. 

The simulations are performed in a vanishing external 
field and we use the fluctuation-dissipation relations to 
evaluate the specific heat and the magnetic susceptibility. 
The averaged value of a quantity X is the result of a time 
average for a sample s followed by an average over n s 
realizations of the disorder so that: 



1 '°s -. i 



(15) 



The total magnetization M and the total energy E are 
defined as 

L L 
M = 22 Ui >i ' E = _ £ ( <J i,j a i+l,] + • 

The quantities studied are the following: 

• The moments of the magnetization per site 



m 



p 



\M\p) 
L 2 p 



m\ = m . 



(17) 



• The susceptibility per site 

AM 2 = (M- (M)f. (18) 
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FIG. 3: Snapshots of equilibrium spin configurations for the 
random system (A = 9) and the pure system (A = 0) in the 
ordered phase (T = 0.8T C ), at the critical point and in the 
paramagnetic phase (T = 2T C ). 
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• The specific heat per site 



ln(L) 



C=±j^, AE 2 = {E-(E)f. (19) 

• The temperature derivative of the logarithm of the 
susceptibility: 

iTr _ JKeam^) AE = E _ {E) _ (20) 

(AM 2 ) 

B. Simulations results 

In Fig.|31equilibrium spin configurations for the system 
with biaxially correlated disorder at A = 9 are compared 
to the corresponding pure system configurations in the 
low-temperature phase, at the critical point and in the 
high-temperature phase for a system with size 116 x 116 
with periodic boundary conditions. In the off-critical sys- 
tems, spin clusters are larger for the disordered system 



FIG. 4: Semilogarithmic plot of the finite-size behavior of the 
specific heat per spin at the critical point. The specific heat 
exponent a is negative and the behavior of C is dominated 
by its regular contribution for large values of L. 



whereas the long-range aspect is roughly the same for 
both systems at the critical temperature. 

The finite-size behaviors at the critical point of the spe- 
cific heat C, the moments of the magnetization density 
m,p, the susceptibility x, and the temperature derivative 
of ln%, \D r , are shown in Figs. 0H7| for three values of the 
disorder amplitude, A = 3, 4 and 9, and for ten sizes 
ranging from L = 8toL=116. 

The specific heat saturates at large system size, which 
indicates that the regular contribution dominates the sin- 
gular one. The singular contribution, behaving as L a l v , 
vanishes at large system size, i.e., the specific heat ex- 
ponent a is negative. This is the reason why we studied 
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.a. 




ln(L) 



FIG. 5: Finite-size scaling of the moments of the magneti- 
zation density at the critical point. The lines correspond to 
linear fits of the data. 
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FIG. 7: Finite-size scaling of ty, temperature derivative of the 
magnetic susceptibility at the critical point. The lines are the 
nonlinear fits of the data as explained in the text. 
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FIG. 8: Linear extrapolation of \jv to infinite disorder 
strength. 



ln(L) 



FIG. 6: Finite-size scaling of the magnetic susceptibility per 
spin at the critical point. The influence of A is quite small 
and the line gives the linear fit of all the data. 



the finite-size behavior of which diverges with the size 
of the system, thus allowing us to estimate the thermal 
exponent 1/v. 

According to finite-size scaling theory, the following 
behaviors are expected for a large critical system with 



size L: 



to p oc L px 
* oc L 1/v , 



(p) 



j/i/ = d- 2x m = 2(1 - x m ) , 
l/v = Vt. (21) 



The log-log plots for (m p ) 1 / p , \, and W display the ex- 
pected linear behavior, although with noticeable devia- 
tions at small size and a slightly A-dependent slope for 
^. Thus the magnetization and susceptibility exponents 
were simply obtained through linear fits, whereas we used 
a nonlinear fit of the data for * 



ln(*) = A + i ln(L) + In (l + BL~") 



(22) 
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TABLE I: Effective critical exponents deduced from the fits 
of the finite-size scaling data in Figs. 00 The last line gives 
the values of the effective correction-to-scaling exponent for 
* defined in Eq. 





A = 3 


A = 4 


A = 9 




0.1312(13) 


0.1292(10) 


0.1294(14) 


T (2) 


0.1283(10) 


0.1279(8) 


0.1309(11) 


T (3) 


0.1259(8) 


0.1264(8) 


0.1310(9) 


T (4) 


0.1238(7) 


0.1249(7) 


0.1305(8) 


j/p 


1.7435(21) 


1.7442(17) 


1.7381(22) 


l/u 


0.5368(9) 


0.5268(10) 


0.5115(8) 




0.737(7) 


0.573(6) 


0.569(7) 



The exponents are collected in Table[|]for the three values 
of the disorder strength. There is no significant variation 
of (3/v and j/u with A. Taking the average of the expo- 
nents obtained for the three values of A leads to 

- = x m = 0.1299(12) , - = 1.7419(20) . (23) 
v v 

There is no clear evidence of multiscaling: The slow decay 

(r>) 

of Xm with increasing p is probably due to a crossover 
effect since it disappears at A = 9 when the disorder is 
sufficiently strong. 

The variation of \ jv with A is significant. As shown 
in Fig. it is linear in 1/A and the extrapolation leads 
to 

- = y t = 0.4988(13) , (24) 
v 

at infinite disorder strength. 

IV. CRITICAL PROFILES 




FIG. 9: Conformal transformation of the half-plane into a 
square. The boundary conditions, free or fixed, are the same 
on both systems. 
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ln(Q 



The analysis of Monte Carlo data for the critical pro- 
files using the tools of conformal theory provides an ef- 
ficient method for the determination of bulk and surface 
critical exponentsiS 4 * 2 ^* 2 ^ 

Using the replica trick for a system with biaxially cor- 
related disorder leads to an effective Hamiltonian with 
long-range interactions between the replicas. Thus we 
do not expect our system to be conformally invariant at 
its new fixed point as it seems to be the case for uncorre- 
cted disorder. 24 Nonetheless, in this section we assume 
that the critical profiles have a leading behavior which is 
the same as for a conformally invariant system. In both 
cases the form of the profiles is constrained by covariance 
under global scale transformations and by the values of 
surface and bulk exponents. We know of at least one ex- 
ample, the random Ising chain in a transverse field, for 
which the critical profiles follow quite accurately the con- 
formal predictions, 27 although the system cannot be con- 
formally invariant since it displays a strongly anisotropic 
scaling behavioriS 2 * 2 ^ 



FIG. 10: Conformal profiles of the magnetization density for 
different system sizes L and disorder strength A. The dashed 
lines indicate the slope —fi/v expected from the finite-size 
scaling results. 



A. Conformal transformation of the densities 

On a semi-infinite critical system with fixed boundary 
conditions, the form of the magnetization profile is fixed, 
up to a constant amplitude, by translational, rotational, 
and global scale invariance. It decays as a power of the 
distance y from the surface^ 



( TO (y))fi XO d = A m y Xm , Xm = P/v , (25) 

where the exponent x m is the scaling dimension of the 
bulk magnetization density. 

Let us consider the Schwarz-Christoffel transformation 
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FIG. 11: Conformal profiles of the singular part of the energy 
density for different system sizes L and disorder strength A. 
The lines correspond to the linear fits of the data. 



of the half-plane z = x + iy such that 
L 



w(z) 



-F{z,k) 



2K(k)' 
where w{z) = u + iv and 



(2Kw\ 
Z = Sn {—) 



F(z,k) = 



dt 

y/(l-t*)(l-&j?j 



(26) 



(27) 



is the elliptic integral of the first kind with modulus k. 
K(k) = F(l,k) is the complete elliptic integral of the 
first kind and sn is the Jacobian elliptic sine. 

This conformal transformation maps the half-plane 
y = ^s(z) > onto the interior of the square with 
-L/2 < u < L/2, < v < L (see Fig. EJ when the 
modulus k verifies the relation K{k)/K(yl — k 2 ) = 1/2 
(see e.g. Ref . l30|) . 

In the new geometry the magnetization density is given 

by 



(™H) fix ed = b(z) = \w'( Z )\- 1 

L 



w'{z) = 



2K(k) y/(l- z 2 )(l~k 2 z 2 ) 
Finally, introducing the reduced variable 

^s[z(w)] 



^{l-z 2 ){l-k 2 z 2 ) 



(28) 



(29) 



TABLE II: Scaling dimension of the energy density x e de- 
duced from the linear fits of the conformal profiles in Fig. 1111 



L 


A = 7 


A = 9 


A = 11 


16 


1.5435(9) 


1.5377(9) 


1.5177(10) 


32 


1.5241(9) 


1.5207(6) 


1.5110(9) 


64 


1.5209(6) 


1.5182(4) 


1.5002(8) 



the magnetization density varies as a power of £ 
m = (m(C)) fixod oc C Xm ■ 



(30) 



The same behavior is obtained for the singular part of 
the energy density, for which 



( e (0)sin E C" 



1/1/ . 



(31) 



In this case one needs to eliminate the regular contri- 
bution to the energy density. The regular part can be 
canceled by taking the difference Se between the energy 
density profiles obtained on a system with either free 
boundary conditions Cf lee (w) or fixed boundary condi- 
tions efi X ed(w) and then 



Se = 



.(C)> 



(efixed(C)) oc C" 



(32) 



The amplitudes of the singular parts have opposite signs 
with the two types of boundary conditions; 26 ! 31 thus in 
the difference the singular contribution is amplified and 
the bulk regular contribution is eliminated. 



B. Monte Carlo simulations 

We work on a square lattice with either fixed or free 
boundary conditions using the Swendsen-Wang cluster 
algorithm 3 ^ which is more efficient than the Wolff algo- 
rithm with fixed boundary conditions. The simulations 
were performed for three lattice sizes (L = 16, 32, 64) 
and three values of the couplings (A = 7, 9, 11). After 
a waiting time r c = 5 x 10 3 MCS, the time averages are 
taken during a time r = 6 x 10 4 MCS. The disorder av- 
erages are taken over a number of samples ranging from 
n s = 2 x 10 4 for the smaller systems to n s = 2 x 10 3 for 
the larger ones. 

For the magnetization density profiles all the surface 
spins are fixed in the +1 state. The behavior of the 
singular part of the energy density is obtained by tak- 
ing the difference of Eq. (|32|) between two systems with 
either free or fixed boundary conditions for the same re- 
alizations of the disorder, which considerably reduces the 
statistical noise. 

The log-log plots for the magnetization and the energy 
densities are shown in Figs. II 01 II II as functions of the 
reduced conformal variable C defined in Eq. (|29|l . 

The magnetization density profiles show strong correc- 
tions to scaling and due to the small values of £ we were 
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TABLE III: Critical exponents deduced from the scaling laws using the finite-size scaling results for /3/v, j/u in Eq. I23H and 
\jv in Eq. (llMI as inputs. 



X m x £ a j3 

0.1294(7) 1.5012(13) -2.01(1) 0.2604(31) 



unable to find a reliable nonlinear fit of our data. Nev- 
ertheless, the dashed lines indicate that the asymptotic 
slopes are in agreement with the value of (3/v obtained 
by finite-size scaling and given in Eq. J22J). 

On the contrary, the corrections are very weak for the 
energy density profiles in Fig. ^] so that the slopes could 
be obtained with a good accuracy through a linear fit of 
the data. The small deviations from linearity occur close 
to the surface for small values of C- Only a few points 
are concerned and their weight in the fit is negligible. 
The values of collected in Table ITT1 for the different 

values of L and A. They are slowly decreasing when A 
or L increases. 



V. DISCUSSION 

Our results are collected in Table ITTT1 The scaling di- 
mension of the magnetization density x rn is directly given 
by (3/v in Eq. 123JI . Alternatively, it can estimated us- 
ing 7/1/ together with the scaling law j/is = d — 2x m , 
which gives x m = 0.1291(10). The value in Tabic InTlcor- 
responds to the centre of the confidence interval which 
is common to the two estimates. It is quite close to the 
value for the pure 2D Ising model x m — 1/8. Thus the 
magnetization in the random and pure systems have al- 
most the same fractal dimensions. This explains why 
both systems in Fig. have roughly the same aspect at 
the critical point. The value of used to calculate 

i] = 2x m + 2 — d = 2x m and 8 = — 1 + d/x m . 

The value of v following from Eq. (|24J) is such that, ac- 
cording to Eq. I|14[) . the disorder is irrelevant at the new 
fixed point, as required for its stability. Furthermore the 
numerical estimate of is, within the error of the calcu- 
lation corresponds to the prediction in Eq. (J2J) with the 
effective correlator exponent, a = 1. The estimate of v 
leads to the value of the specific heat exponent through 
the Josephson scaling law a = 2 — dv. The correlation 
length diverges more strongly in the random system with 
v close to 2 than in the pure system with v = 1, this 
might explain the difference in the cluster sizes for the 
off-critical systems in Fig. [21 

The scaling dimension of the energy density x e is ob- 
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3.492(13) 14.46(8) 2.005(5) 0.2588(14) 



tained using the scaling law x e = d— \ jv. The values of 
x e deduced from the conformal profiles, collected in Ta- 
ble[H] are close to this estimate for large sizes and strong 
enough disorder. 

To summarize we have studied the critical behavior of 
the square lattice Ising model with biaxially correlated 
disorder by large scale MC simulations. The different 
estimates of the critical exponents give a consistent pic- 
ture, in which v seems to be the same as for isotropic 
Gaussian correlated disorder, whereas x m is not very far 
from the pure and uncorrelated random systems values. 

One important technical simplification of the model is 
its self-duality property, which is used to locate exactly 
the critical point. The same type of parametrization, as 
given in Eq. JSJ) can be used to preserve self-duality in 
the q-state Potts model, for any value of q~ An inves- 
tigation of these models could provide further informa- 
tion about the nature of the phase transition in disor- 
dered systems. For example, q = 1 is equivalent to the 
percolation problem^ whereas the phase transition in 
the pure system is of first order for q > 4. 35 In the lat- 
ter problem, due to uncorrelated disorder, the transition 
is softened into a second-order one^ with a correlation 
length exponent v s h or t « 1 — Consequently, according to 
the extended Harris criterion in Eq. Q), the phase tran- 
sition is expected to belong to a new and g-dependent 
universality class. If the limiting behavior for large q can 
be predicted, as for uncorrelated disorder— remains the 
subject of a separate study. 
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